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Abstract. In this paper, we discuss simulations of 
MHD wave generation and propagation through a three- 
dimensional open magnetic flux tube in the lower solar at- 
mosphere. By using self-similar analytical solutions for 
modelling the magnetic field in Cartesian coordinate system, 
we have constructed a 3-D magnetohydrostatic configuration 
which is used as the initial condition for non-linear MHD 
wave simulations. For a driver we have implemented a high- 
frequency vortex-type motion at the footpoint region of the 
open magnetic flux tube. It is found that the implemented 
swirly source is able to excite different types of wave modes, 
i.e. sausage, kink and torsional Alfvén modes. Analysing 
these waves by magneto-seismology tools could provide in- 
sight into the magnetic structure of the lower solar atmo- 
sphere. 


Keywords. Solar physics, astrophysics, and astronomy 
(Photosphere and chromosphere) 


1 Introduction 


There is observational evidence of a variety of magnetic field 
structures, i.e. open magnetic flux tubes, solar coronal loops, 
etc. in the higher layers of the solar atmosphere. These 
configurations, while in magnetohydrostatic equilibrium and 
thus relatively long-lived, are subject to wave motions, if 
a perturbation acts on them. The studies of recently re- 
ported observed ubiquitous wave motions in a range of mag- 
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netic field configurations in the photosphere, chromosphere 
and corona are of particular interest, since through different 
mechanisms of wave damping they might be responsible for 
the heating of the solar plasma (see e.g. Banerjee et al., 2007; 
Taroyan and Erdélyi, 2009). Furthermore, one can exploit 
these waves by the technique of solar magneto-seismology to 
probe the fine structure of the Sun’s magnetised and highly 
dynamic atmosphere (see e.g. Goossens et al., 2002; Ar- 
regui et al., 2007; Erdélyi and Fedun, 2007; Goossens et al., 
2008; Verth, 2008; Andries et al., 2005, 2009; Ruderman and 
Erdélyi, 2010). A popular excitation mechanism of waves in 
magnetic chromospheric and coronal structures is the leak- 
age of oscillatory modes from the inner layers of the solar 
photosphere along magnetic regions, because the magnetic 
field favourably shifts the acoustic cut-off frequency. 

Periodic motions at the footpoint regions of magnetic flux 
tubes can generate waves that may supply energy to the upper 
part of the solar atmosphere. For example, granular buffet- 
ing motion could be responsible for the excitation of kink 
(transverse) waves. Turbulence in the convection zone near 
the surface of the Sun excites solar p-modes (Goldreich and 
Keeley, 1977). These acoustic waves are transmitted to the 
photospheric region of the solar atmosphere and are able to 
drive longitudinal magneto-acoustic waves in magnetic flux 
tubes. These motions have large spatial scales. A variety 
of one- and multi-dimensional MHD numerical simulations 
of wave generation in a solar magnetic flux tube by various 
periodic vertical or horizontal sources have already been per- 
formed (see for example Bogdan et al., 2003; Hasan et al., 
2005; Vigeesh et al., 2009; Ofman, 2009; Fedun et al., 2009; 
Felipe et al., 2010; Fedun et al., 2011). Here we push these 
frontiers further. 
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Fig. 1. Initial magnetic flux tube configuration. Selected colored 
lines correspond to the magnetic field lines. Also, color indicates 
the value of the magnetic field strength at different regions of the 
flux tube. The iso-contours of the initially constant magnetic field 
are overplotted at the top of the simulation box. Note, only central 
part (i.e. x, y = [0.8, 1.2] Mm) of the full computational domain (i.e. 
x,y =[0,2] Mm) is shown. 


Modern state-of-the-art ground- and space-based solar ob- 
servational instruments e.g. Swedish Solar Telescope (SST), 
DST/ROSA, Hinode, SDO etc., carry out the highest spatio- 
temporal resolution of dynamical processes in the lower and 
upper regions of the solar atmosphere (see e.g. Bonet et al., 
2008; Wedemeyer-B6hm and Rouppe van der Voort, 2009; 
Wedemeyer-Bohm, 2010). 

Recently, Bonet et al. (2008) have discovered localised 
vortex-type motions created at downdrafts where the plasma 
returns to the solar interior after cooling down. The reported 
swirls have diameters 0.2-1.5 Mm and lifetimes of 3-7 min. 
They have found numerous examples of convectively driven 
vortex flows with clockwise and counterclockwise rotation. 
Wedemeyer-Bo6hm and Rouppe van der Voort (2009) have 
analysed the time series of spectral scans through the Ca II 
854.2 nm spectral line with the CRISP instrument mounted 
at the Swedish Solar Telescope and found rotations up in 
the chromosphere. They interpret these swirly motions and 
the associated bright point motions as a direct indication of 
upper-atmospheric magnetic field twisting and braiding as a 
result of convective buffeting of magnetic footpoints. We 
show that such observed rotational motions could be a natu- 
ral driver for not just Alfvén but a range of MHD waves, i.e., 
slow (SMAW) and fast (FMAW) magnetoacoustic waves. 

Previously, the torsional Alfvén wave has been studied 
theoretically by (e.g. Hollweg, 1981; Hollweg et al., 1982). 
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In these papers the authors have analysed the linearised prop- 
agation of axisymmetric twists on axisymmetric vertical flux 
tubes. Open and closed configurations of the magnetic flux 
tubes which may model e.g. the coronal hole and active re- 
gion loops have been studied. It was found that torsional 
Alfvén waves may produce enough energy flux to heat the 
solar corona (see also Antolin et al., 2008). 

Verth et al. (2010) have shown that observation of the 
eigenmodes of torsional Alfvén waves can provide tempera- 
ture diagnostics of both the internal and surrounding plasma, 
i.e. these waves are the ideal magneto-seismological tool for 
probing radial plasma inhomogeneity in solar waveguides. 

Due to the incompressible nature of the Alfvén wave the 
detection of this type of wave MHD mode is rather challeng- 
ing and is difficult in the solar atmosphere. Only recently 
Jess et al. (2009), by studying of a chromospheric magnetic 
bright point (MBP) group and analysing periodic Hæ spec- 
tral line broadening with no intensity variations, have clearly 
shown the existence of this wave mode in the Sun. 

Numerically, vortex-type motion in the radiative MHD 
simulations of magnetoconvection has recently been found 
by Shelyag et al. (2011a,b). These authors have shown a 
direct connection between magnetic vortices and rotary mo- 
tions of photospheric bright points, and suggested that there 
may be a connection between the MBP rotation and small- 
scale swirly motions observed higher in the atmosphere. 

Based on these theoretical, observational and numerical 
results, in this short paper, we focus on waves emanating 
from a spatially localised vortex source at the footpoint re- 
gion of an open magnetic flux tube. Due to the swirly mo- 
tion it is anticipated that such drivers are the natural exciters 
of torsional Alfvén waves. Analysing the evolution of the 
horizontal cross section of a simulated magnetic flux tube 
we are able to establish clear evidence of propagating MHD 
kink and sausage modes as well. 


2 Numerical model and simulations 


We use the 3-D version of Sheffield Advanced Code (SAC) 
to perform our simulations (for more details about code see 
Shelyag et al., 2008). The computational domain has a phys- 
ical size of 2 x 2 Mm in the horizontal directions x and y, 
and in the vertical direction z the domain extends from the 
lower photosphere through the chromosphere up to the tran- 
sition region i.e. 1.6 Mm in height. Numerically the domain 
is resolved by 100 x 100 x 196 grid cells in the x-, y- and z- 
directions, respectively. The boundaries of the computational 
box are set to be open. The standard VAL IIC (Vernazza 
et al., 1981) model profiles of temperature and density have 
been used as a realistic stratified hydrostatic background. 
The self-similar approach is used to construct the non- 
potential 3-D magnetic field configuration, i.e., open mag- 
netic flux tube. Schlüter and Temesvary (1958) were 
the first who suggested to use self-similar magnetic field 
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Fig. 2. Rendering of the magnetic field (by, bg, bz) and velocity (Vr, Vg, Vz) perturbation components generated by the vortex-type periodic 
driver (see Eq. 2) at the footpoint of the magnetic flux tube. Two different times are shown. The top and bottom sets of six horizontal slices 
at hight h =0.12 Mm correspond to simulation times t = 1.5 s and t = 120.02 s, respectively. 


constructions for description of various solar magnetic phe- B, =— * d Boz (z) -G(f), 

nomena. Currently, magnetic fields of this type are widely 2 dz 

used in numerical simulations in helioseismology and coro- AE y d Boz (z) Gf), (1) 
nal physics (e.g. Gordovskyy and Jain, 2007; Cameron et al., 2 dz 


Boze) G(f), 


2008; Shelyag et al., 2009; Fedun et al., 2011). In our par- B, 
ticular three-dimensional case, due to the self-similarity, the 

open magnetic flux tube is obtained analytically from the fol- 

lowing set of equations: 
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Fig. 3. Snapshots showing the iso-surfaces of the horizontal velocity component Vx during numerical simulation. From left to the right 
the snapshots are at t = 7.2, 312s at the centre of the numerical domain we have overplotted the 2-D vertical slice of the background Bz 
component of the magnetic field. Note, rendering of the 3-D numerical data has been constructed using the VAPOR visualisation package. 
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Fig. 4. The time-distance diagrams of the radial (V,) and the vertical (Vz) components of the velocity at x = 0.98 Mm and y = 0.98 Mm, i.e., 


near the axis of the magnetic flux tube. 


and 
x? + y? 
2 
ro 


J= Boz (2), (2) 
where Bọoz(z) is the vertical z component of the background 
magnetic field along the symmetry axis towards the top 
boundary of the model, rọ is a radial scaling and G is an 
arbitrarily chosen function which describes the profile in the 
radial direction of the vertical magnetic field component. In 


the present case this function is chosen to be exponential, i.e. 


G(f)=A-exp(—f), (3) 


where A is an arbitrary constant. To implement this magnetic 
field configuration into the background hydrostatic equilib- 
rium, we have used the same principles as Fedun et al. (2011) 
for a two-dimensional case. The constructed 3-D magnetic 
flux tube is then visualised by selected magnetic field lines as 
demonstrated in Fig. 1. The different colors along the field 
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lines correspond to magnetic field strength. The magnetic 
field strength at the footpoint region is Bo = 1000 G and re- 
duces with height in the mid-chromosphere to a few G (in 
the horizontal direction the magnetic field has a maximum 
absolute value at the centre and minimum value at edges of 
the flux tube). The footpoint flux tube radius is 100 km. 

In the Introduction we have noted that vortex-type motion 
at the photospheric region is widely observed. Swirly mo- 
tions and associated MBP motions are a direct indication 
of upper-atmospheric magnetic field twisting and braiding 
as a result of convective buffeting of magnetic footpoints. 
Based on these observational results we have implemented 
here a vortex-type high-frequency periodic source located at 
the footpoint of the magnetic flux tube as an initial driver in 
our numerical simulations. The Vx and Vy components of the 
velocity perturbation have a Gaussian spatial distribution in 
the x-, y- and z-directions: 
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Fig. 5. A time series of horizontal cuts at height h = 0.75 Mm of the total magnetic field. The strength of the magnetic field is shown as 
iso-contours. The difference in time between the snapshots is approximately equal to 30 s, i.e., the period of footpoint driver. Note the clear 
evidence of the sausage and kink oscillations which propagate along the magnetic flux tube. 
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where Ag is the amplitude of the initial perturbation; Ar and 
Az are the half width of the Gaussian spatial profiles of the 
driver in the radial and vertical directions, respectively; T is 
the period of the driver, r? = x? + y? is the radial distance. In 
our numerical simulations Ar = 0.1 Mm and Az =0.01 Mm 
in order to mimic observed drivers. 

In Fig. 2 the initial state and disturbances generated by a 
vortex driver at a later elapsed time (t = 120s) are shown as 
a two-dimensional horizontal slice at height h = 0.12 Mm, 
i.e. just above the region where the driver is located in the 
computational domain. We have plotted a color rendering 
of the perturbed components of the magnetic field (br, bg, 
bz) and velocities (V;, Vg, Vz) in cylindrical coordinate sys- 
tem. Disturbances excited by the implemented vortex driver 
in the computational domain are similar to those obtained 
by Shelyag et al. (201 1a,b) in direct radiative MHD simula- 
tions of magneto-convection and observed vortex motions by 
Bonet et al. (2008) in photospheric G-band bright points. 

The horizontal component of velocity Vx is shown in Fig. 3 
as isosurfaces. Illustrative snapshots are captured at times 
t =7.2 and 312s, respectively. At the centre of each snap- 
shot we have overplotted the vertical cut of the background 
magnetic field component B, which shows the internal mag- 
netic structure of the constructed magnetic flux tube. The 
high frequency (T = 30s) vortex driver is located at the foot- 
point region of the magnetic flux tube (see the first upper left 
panel of Fig. 3). 

After ten periods (see the right panel of Fig. 3) we ob- 
serve a 3-D pattern of the SMAW and FMAW distribution in 
the computational domain. The SMAWs are located near the 
axis of the magnetic flux tube and propagate vertically up- 
ward, while the FMAWs propagate obliquely. It can be seen 
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that SMAWs, generated by the driver, are weak and present 
only up to the height A = 0.7 Mm, i.e., SMAWs do not supply 
considerable energy to the upper part of the chromosphere. 
Note, in a 3-D geometry we can now clearly see the shape of 
generated waves. SMAWs have a maximum of amplitude at 
the axis of the magnetic flux tube, equal to 800 ms7!. 

Time-distance diagrams of the radial (V,) and vertical (Vz) 
components of the velocity rendered approximately at the 
axis of the magnetic flux tube are shown in Fig. 4. As 
we have noted previously, two types of waves with differ- 
ent phase speeds can be resolved. The SMAW and FMAW 
propagate with phase speeds Vph = 3.8 km s7! and Vph = 
8.2 km s71, respectively. 

The dark horizontal regions on the time-distance diagram 
of the radial component (V,) (see left panel of Fig. 4) occur 
due to wave interference (both SMAW and FMAW) in the 
radial direction. 

In Fig. 5 we have shown iso-contours of the total mag- 
netic field at h = 0.75 Mm for simulation times t ~ 142, 171, 
199 and 227. It is clearly seen that the area inside, e.g., 
the iso-contour at 32 Gauss, decreases and increases period- 
ically. This is clear evidence of a sausage mode excited by 
a vortex driver. Could this be observed? Our prediction of 
generated sausage waves has recently been confirmed obser- 
vationally by Morton et al. (2011). Furthermore, this driver 
also excites the kink mode as is shown by transverse motion 
of the magnetic flux tube centre (see Fig. 5). We appreciate 
that is difficult to identify the changes of position and areas 
of iso-contours in static images. To help readers, we also 
provided movies in downloadable electronic form, that can 
be found at http://swat.group.shef.ac.uk/simulations.html. 


3 Conclusions 
The purpose of this short report is to investigate numeri- 


cally the generation of MHD waves in an open magnetic flux 
tube by a vortex-type photospheric driver in a very realistic 
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three-dimensional geometry. We have shown that a such 
driver can excite both SMAW and FMAW. The torsional 
Alfvén wave is also present and the interesting properties 
of this mode will be a focus of our future study. There is 
simply no space available here for a deep analysis. We have 
not investigated the possible heating contribution of the gen- 
erated MHD modes. This requires much more detailed nu- 
merical and analytical analysis. For more efficient heating 
due to MHD waves, these waves must transfer their energy 
to smaller length scales due to various types of instability 
and dissipation mechanisms (see Fedun et al., 2004; Copil 
et al., 2010). In order to resolve the dissipation numerically 
very powerful high performance computational equipment is 
needed and we hope to make the first steps in this direction 
in the near future. 
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